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Abstract 

We  consider  the  use  of  a  harmonic  oscillator  model  to 
investigate  the  tunnelling  of  atoms.  Ions,  and  molecular  groups 
between  positions  of  equilibrium)  for  systems  in  condensed  phases. 
The  use  of  the  oscillator  model  is  warranted  in  many  instances 
for  which  one  has  reasonably  accurate  information  about  the 
vibrational  frequencies  of  atoms  or  molecular  groups  together  with 
information  about  the  distance  over  which  the  species  transfers 
Thus,  the  harmonic  oscillator  model  can  be  used  in  the  absence 
of  more  complete  information  to  make  estimates  of  the  rates  of 
various  kinds  of  heavy-centre  transfers  in  solution.  In  this 
paper  we  present  a  general  treatment  for  the  determination  of 
the  tunnel-factors  which  are  involved  in  the  transfer.  The 
work  constitutes  an  extention  of  Che  analysis  of  Wall  and  Glockler 
who  used  a  symmetric  double  potential  well  in  order  to  investigate 
the  resonance  doubling  of  the  vibrational  spectrum  of  anmonia. 

We  consider  the  general  problem  of  the  tunnelling  between  any 
two  vibrational  states  of  the  separate  potential  wells.  As 
an  example,  we  apply  the  analysis  to  the  consideration  of  a  simple 
ion- transfer  between  two  centres  of  solvation. 
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In  t  ruduet ion 

The  transfer  of  an  atom,  ion  or  molecular  group  from  one 
location  to  another  in  a  condensed  phase  was  considered  in  Part  l.^ 
There,  methods  were  developed  which  make  use  of  accurate 
pair-wise  potential  energy  functions  to  determine  the  magnitudes 
of  the  matrix  elements  which  govern  the  transfer.  For  solution 
chemistry  it  is  clear  that  it  may  not  generally  be  possible  to  use 
accurate  potential  energy  functions,  as  few  accurate  functions  have 
been  determined.  Thus,  it  is  difficult  to  construct  rate  constants 
and  to  establish  links  between  structural  and  kinetic  factors 
of  the  kind  considered^  without  corroborating  experimental  data 
of  different  kinds. 

In  this  paper  we  consider  the  use  of  the  relatively  simple 
double  harmonic  potential  well  in  order  to  describe  the  parametric 
interaction  between  a  transferable  species  at  one  location  and 
its  eventual  location  at  a  physically  removed  side.  We  show  that 
it  is  possible  to  make  an  estimate  of  the  size  of  the  matrix  elements 
which  couple  the  initial  and  final  Btates  of  the  system.  These 
estimates  can  be  made  with  only  the  frequencies  of  vibration  of  the 
transfer-species  in  its  initial  and  final  states  and  the  distance 

over  which  the  species  tunnels.  As  the  analysis  of  Wall  and 
2 

Glockler  on  the  inversion  tunnelling  of  ammonia  showed,  the  height 
of  the  energy  barrier  between  the  initial  and  final  states  can  be 
aa  much  as  507,  in  excess  of  the  true  value.  Yet,  with  the  use 
of  simple  perturbation  theory,  Wall  and  Glockler' s  analysis  is  sur¬ 
prisingly  accurate.  The  height  of  the  barrier  for  the  amnonia 
inversion  has  been  accurately  estimated  by  Swalen  and  Ibers.^ 


-  3- 

In  this  paper  the  Wal  1  -Glockler  analvsis*  i  *.  «xt*nlt*d  .vui 
applied  to  the  general,  asymmetric,  harmonic  double  potential 
well.  The  results  of  the  generalization  are  applied  to  the 
determine  .ion  of  the  rate  constant  for  the  tunnel  1  in?  of  a  heaw- 
centre  from  one  location  co  another.  We  present  a  general  analysis 
which  parallels  that  of  Part  I.  The  specific  difference  between 
this  work  and  that  of  Part  1  is  the  use  of  the  double  potential  well 
instead  of  considering  the  specific  pair-wise  potential  energv 
functions  which  operate  between  the  transfer- species  and  its  surround¬ 
ing  solvent  or  ligand. 

In  the  next  section,  we  sunxnarize  the  Born- 

4 

Oppenheitner-Holstein  separation  which  was  developed  in  Part  I 
Following  that  stannary.  we  present  specific  forms  for  the  matrix 
elements  which  arise  from  the  considerat ion  of  the  double  potential 
well.  Finally,  the  results  are  applied  to  a  model  type  of  ion 
transfer . 


penheimer-Hol stein  Analysis:  Stannary 


The  Hamiltonian  operator  for  the  entire  reactive  system 
is  expressed  as 


in  which  Ht  is  that  part  of  the  total  operator  which  applies  to 
the  transfer- species  and  Henv  Is  the  remainder.  The  operator  Ht 
is  expressed  further  as 
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in  which  T  is  the  momentum  operator  and  V  is  the  potential. 

These  operators  usually  can  be  specified  in  one  dimension  for  the 
motion  of  the  transfer- species  along  the  transfer-axis.^ 

The  complete  eigenvalue  problem  is  written  as 


in  which  the  basis  functions  ^  are  harmonic  oscillator  functions 
which  are  referred  to  coordinate  origins  at  the  locations  i 
(i  •  a  or  b,  the  labels  for  the  Initial  and  final  locations  of 
the  transfer-species).  The  elements  X^  of  the  vector  X  are 
in  essence  the  vibrational  wavefunctions  for  the  environment . 

A  particular  element  of  the  vector  X  is  determined  by  the 


H.  <  xi  -  l  41X., 

in  ,in  i>  ^  j6 


In  which  the  L-matrlx  elements  are  defined  by 


The  S-  ^  are  elements  of  the  inverse  of  the  overlap  matrix 
The  diagonal  matrix  element  ^  is  given  by 


in  which  E^  is  the  eigenvalue  of  the  harmonic  oscillator  Hamiltonian 
operator  for  the  transfer-species  at  the  location  i.  Finally, 
the  potential  energy  operator  ,  which  enters  the  matrix  elements 
Vkc  j s  in  e<Jn  (6),  is  defined  as 


Vj  -  V  -  fckzj 


for  which  2  is  defined  to  be  the  transfer-axis,  and  kj  is  the 
harmonic  oscillator  force  constant  which  is  defined  with  reference 
to  the  coordinates  at  the  origin  j.  This  form  for  the  potential 
is  used  in  order  to  isolate  the  local  harmonic  oscillator 
Hamiltonian  operators.1 

In  the  next  section,  the  matrix  elements  V.^  .  .  which  contribute 

Kt  ,  j  s 

to  the  eqn  (6),  are  determined  for  the  general  asymmetric 

double  harmonic  potential  well. 


Matrix  Elements:  evaluation 


The  general  form  of  the  asymmetric  double  potential  well  is 
illustrated  in  Figure  1.  The  actual  potential  is  specified  as 
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is  used  in 
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-  *Vb  Oi21-  (9) 

in  which  -V0  is  a  constant  displacement  in  energy  of  the  left  hand 
well  relative  to  the  right.  In  terms  of  a  Taylor  series  expansion, 
-V(,  represents  the  difference  Vb(z£)  "  Va^2a^ ’  t*ie  difference 
between  the  separate  zeroth  order  terms  of  the  expansion;  note.  V0 
can  be  negative,  in  which  case  the  right  hand  well  is  lower  in 
energy  than  the  left. 

The  intersection  of  the  two  wells  is  defined  to  be  the  origin 
of  the  system  of  coordinates.  The  distances  from  the  origin  to  the 
minima  of  the  two  wells  are  ia  <to  the  left)  and  (to  the  right). 
If  is  known,  for  example,  then 

*»  ’  l(kb/k«),b!  +  2V»/ka,%'  <10) 


+  %kbj  dz  H.m(z,).nUb>  -  *kbf  iz  *i*a“.>*n(V 
0 

-  •  kbIb  '  2v0»W‘n(*b>  (13) 

The  specific  form  is  used  here  for  the  simple  reason  that  there 
are  only  two  centres,  a  and  b,  and  the  exchange  of  indices  converts 
one  to  the  other. 

The  integrals  in  eqn  (13)  are  evaluated  with  the  use  of 
the  one- dimensional  harmonic  oscillator  basis  functions 


With  R  -  i  +  i.  ,  one  finds 
a  b 

*.  *  2  K/k.  +  2<2  -  Vk.>v°/k.R)']  ■  <U) 

However,  if  k.  •  k.  .  then 

a  d 

ta  -  *R  +  V0/kR.  (12) 

For  the  evaluation  of  the  integrals,  the  distances  and  R 

are  specified  Individually.  The  relationships  between  these  } 

quantities  are  as  given  in  eqn  (10)  to  (12). 

The  modified  potential  energy  operator  (j«a,b),  eqn  (8), 


♦n(z)  •  (aV2n«nl)*Hn(ar)exp(-%aJz2) 


(K) 


in  which  Hn(az)  is  the  Hermite  polynomial.  In  the  treatment  of 
the  general  form  of  matrix  element,  it  is  necessary  to  refer  the 
product  of  the  wavefunctions  to  a  single  centre.  This  is 
accomplished  as  follows.  The  exponential  part  combines  as 


exp{-4(a2z2  +  b2z£) > 


exp  -%(a2+b2)y2 


a?b? 

2 (a2+b2) 


R7 


(15) 


with  y  defined  as 


y  •  x  +  d. 


(16) 
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a-  ♦  b- 


For  the  sing ' e-cent re  matrix  element,  R  -  0,  end  eqn  (15)  reduces 
to  exp<  -a-'z-  ) 


The  matrix  element,  eqn  (13),  can  written  as 


.  J  V  1  bn  >  •  4  I 


4S00(-l)”m’|  dy((ka  -  kb)y2 


-9- 

with  p  -  0,1,2,  must  be  evaluated.  The  addition  theorem  for  the 
Hermite  polynomials  is"5 


V.(y-.R)]  ■  (22> 

where  (™)  is  the  binomial  coefficient.  Further,  the  general  form 
of  the  Hermite  polynomial  is^ 


HJx)  -  lm(2)(-l)*(2x)“-2* 
s“0 


ml 

rm-2s)ITT 


(23) 


-  2<“A  +  W*  +  kA  -  SA  -  2V«J 

■H.Dfa(y-Rb,)Hn(b^+Ra)i*xi,(-!i<*2+b?)yJ) 


(18) 


in  wnich  S0o  is  the  ground  state  overlap: 


a2b2 

— -  exp  -  — - R2 

a‘+b2j  2(a2+b2) 


(19) 


a  a2+b2 

Rh  -  — R.  (20) 

b  .!+b' 

It  is  apparent  that  the  general  integral 

f  dy  yPHn,(«(y-Rb))Hn(b(y+Ra))exp(-%(«2+b2)y2)  (21) 

'  -d 


the  notation  [m/2]  indicates  that  the  integer  of  m/2  is  to  be 
taken. 

With  the  use  of  the  expansions,  the  matrix  element  can  be 
expressed  as 


<xm|Vb|bn>  -  -5iTj(.x.lnl)-'st,t,  l  l  (£>(?)(-l)n'k 
2  k-o  1-0  1 


’Hn.-k<^Rb)Hn-.<-,2bV22(k+')[-^r7)%<k+,)k 


lk(2J  l‘(2>  (-l)*+t 

i0  tio  nc-2s)l(i-2t)!s!t 


a2+b2 

af+b2 


4/2a2b2 


<k.  -  V 


_i_  -  2(kA  ♦  Wl-i- 


lh±*tl  -s-t.D2  +  (kaR^  -  kbR2  -  2V0)r  fk*g*1-  «a-t.D2 


(24) 


where  D2  •  4(*2+b?)d2  end  the  Incomplete  gamma  function  is3 


r(a,x)  -  I  dt  c*"1*"' 
•>x 


(25) 


A 

4* 


i 
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The  single-centre  matrix  element  can  be  obtained  in  principle 
from  the  tvo-centre  form,  but  it  is  actually  simpler  to  consider 
the  specific  form 


-2t. - (k+ 1 -  2  (s+t )  - 1 ) !J_ 


'Va'an'  '  11  (2^~m.nl)1»l-dZ  '  Vb  -  2Vc>Hm(«.) 


(26) 


and 


k+t -2 (s+t)  *  even. 


%  |m/2|  [n/2 J 

l  l  <-l)S+t  4  8  t 
s-O  t-0 


The  same  type  of  analysis,  as  applied  to  the  two-centre  form,  yields 
*  am  |  Va  |  am  - 

x(tr,-2sH“n-5tVleltl[(1  ‘  VVl^TT2  "**t  •e»i] 

-  2(kb/k<1)«Rr(Ht^  -s-t,q»)  -  ( (k,,/^)  («R)  2 


+  2V0/Kua  -.-c.q*]  (27) 


for  which  q#  is 


qb  is  defined  similarly. 

Finally,  the  general  overlap  Integral  la  written  aa 


Son  '  j0  M-k<'J»VV.<-7bR. 


m-k 


|k(Jl  <‘<2>  a... 

,£0  ti0  TV-'7i)TT.^tTr.-|tl(-U  (J,7> 


kill 


■+t,,  «,k  +i-2  (a+t) 


Model  System  Calculation.  Asynaetric  ground-to-ground  state  transfer 

The  simplest  direct  application  of  the  methods  which  have 
been  developed  above  is  to  the  consideration  of  the  ground-to- 
ground  state  transfer  of  a  heavy-centre  from  one  harmonic  well  tc 
another  The  wells  need  not  be  symmetric  with  respect  to  their 
widths  nor  with  respect  to  their  relative  locations  on  the 
energy-axis.  According  to  the  transition  state  theory,  the 
convolution  of  environmental  energies  ensures  the  conservation  of 
total  energy  for  the  transfer;  the  energy  deficit  for  the  reactive 
subsystem  is  made  up  from  contributions  from  the  environment 

In  this  section  attention  focuses  entirely  upon  the  problem  of 
determining  the  sire  of  the  matrix  elements  b0  in  the  diabatlc 
limit.  For  the  ground  state  transfer,  we  need  specifically  to 
consider  the  element 

Ls0,b0  “  *Va0,b0  ‘  S°#Vb0tb0^1  *  S°o)‘  (29) 

The  matrix  element  b0  is 


1 

j. 


•  W“hso»1rT 


r(l  +  k  /k.  )rO/2.D-') 


:f - ^ - 1  bt(k  /k.)R.  +  R  ir(i.D-)  +  )  (bR  ) : 

U+a-/biJ  a  6  -  a  * 


(ka/kb)(bRb)J-  2V0/Rwb)r<t,D-) 


l+a^/bM  1+a- /b2  I 


Vb0.b0  *  'bO|VblbO» 


i  H-b  (ka/kb  -  l)r<3/2,q£>  -  2(k4/kb)bRr(l,q^) 


+  (ka/kb)(bR)J  -  2V0/JUb  r(4,q=) 


The  vibrational  overlap  element  for  the  ground  states  is  given  by 
eqn  (  19)  • 

The  general  formulae  for  the  matrix  elements  can  be  simplified 
in  many  instances  with  the  use  of  various  limiting  forms.  The 
limits  generally  involve  the  consideration  of  large  or  small 
values  of  the  arguments  of  the  incomplete  gamna  functions. 

The  typical  values  of  ta,  and  the  several  R’s,  together  with  the 
coefficients  a  and  b.  lead  to  terms,  c.g->  qa .  D.  eta.,  which  are 
large  or  small  in  a  limiting  sense. 


For  integer  n,  the  incomplete  gamna  function  :(n,x->  is 


r(n,x2)  *  (n-l)1e"x  l  x^s/sl 
s-0 


while  for  half-integer  arguments,  n+fc ,  one  has 


r(n+4.,2)  .  iinziUlU  erfc<*>  ♦  X 


in  which  erfc(x)  is  the  complement  of  the  error  integral 


.  ±  r  dt 

a  Jx 


1  -  erf(x). 


In  the  limit  as  x  vanishes. 


lim  r (a.x2)  -  r(a)  . 
x-0 


Conversely,  as  x  becomes  large  (x  >  3).  erfc(x)  -  0.  and 

/ | v i |  _jH*l  ,S+l„2x+l 
lto  rlnH.t1)  -  (  -  o' 


The  application  of  these  limiting  forms  to  the  matrix  element* 
results  in  a  considerable  simplification  In  order  to  apply  these 
limiting  forms,  however,  it  is  necessary  to  consider  the  physics, 
distinctions  within  the  system  which  will  govern  the  types  of 
limitB  which  one  may  expect  to  use. 


1 
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A  nur.ln  r  of  possible  limits  can  be  c^'sidiri'd,  these  limits 
aie  i  f  ied  through  the.  consideration  of  differences  in  the 

reduced  masse:,  in  the  initial  and  final  states,  differences  in 
the  lengths  t  and  and  differences  in  the  force  constants. 

For  the  consideration  of  the  transfer  of  a  simple  solvated  ionic 
species  (for  example,  a  sodium  cation)  from  one  cage  of  solvent 
to  another  in  a  homogeneous,  single  solvent  system,  it  is  clear 
that  (effectively)  all  quantities  in  the  final  state  are  the  same 
as  in  the  initial  state  The  double  harmonic  well  representation 
is  symmetric.  Thus,  clearly,  the  quantity  D  •  0.  On  the  other 
hand,  it  is  possible  to  use  the  harmonic  oscillator  representation 
developed  here  to  consider  the  matrix  element  for  the  cransfer 
of  a  simple  cation  (again,  sodium)  from  one  solvent  system  to 
another.  In  this  case,  as  we  now  show,  we  can  still  obtain  a 
relatively  simple  limiting  form  for  D  and  the  other  quantities. 

We  consider  the  hypothetical  transfer  of  a  simple  cation 
across  the  interface  of  two  polar,  but  insni scible  liquids,  for 
example,  across  the  water-propylene  carbonate  interface.  Although 
in  any  interface  there  do  not  occur  discontinuous  changes  in 
composition,  there  may  be  auff icier t lv  different  structures  of 
solvation  at  some  point  to  ensure  that  the  ion^transfer  between  these 
structures  is  rate-determining.  We  model  such  a  situation. 

The  purpose  in  using  the  harmonic  oscillator  well  model,  as 
stated,  is  to  be  able  to  make  reasonable  estimates  of  the  transfer 
matrix  elements  with  the  use  of  available  data.  The  data  which 
are  needed  in  order  to  consider  the  interfacial  ion-transfer  are 
simply  the  fundamental  frequencies  for  the  oscillation  of  the 
Ion  In  its  initial  and  final  cages  of  solvent  and  the  distance 


over  which  tin-  transfei  t  at  i  •  plait 

Studies  of  the  far  -  inf  rated  spectra  >' f  s'lvated  cation 
of  Croup  I.  sodium.  -  •  ■  .  indicate  a  relatively  small 

solvent- dependence .  '  that  is,  the  change  from  n:\e  solvent  ►  anetbet 
resells  in  o  shift  of  the  spectral  maximum  '  rrlarivclv  small 

amount  The  changes  in  the  reduced  mass  ’  rsir.p  from  on c 
solvent  to  another  are  small  a?  well  Th'  v  a:;  nr  t  be 

reasonable,  as  a  first  approximat  ion .  to  r  t!'«  differences 

in  masf^  In  considering  the  transfer  matr.  's  Moreover 

if  one  can  ignore  the  differences  in  mass.  ’•  .  as  observed.' 

the  spectra  in  several  solvents  are  essentially  the  same,  then 
it  Is  also  reasonable  to  assume  that  the  fore,  constants  are 
effectively  similar  for  the  initial  and  final  states  of  the  ion 
transfer.  With  the  ust  of  these  two  (apparently  extreme'  assumption 
it  is  easy  to  show  that  D  reduces  to 


D  -  V0/2h*aR  (37) 

in  which  w  is  now  the  cotiwon  frequency  and  P.  is  still  the  total 
distance  for  the  transfer  V0  characterizes  the  relative  displace¬ 
ment  of  one  potential  energy  function  (from  which  the  second 
derivative  yields  the  force  constant)  with  respect  to  another. 

Thus,  this  quantity  must  necessarily  reflect  the  total  interaction 
of  the  ion  with  its  surroundings.  As  such,  V0  is  effectively 
the  change  In  free  energy  for  the  ion  in  one  cage  of  solvent 
compared  to  the  other.  ^ 

If  we  assume ,  for  sodium,  that  u.  ~  200  cm  ,  and  R  O.b  A 
we  find 


it  i.*.  apparent  that  large  differences  in  Che  free  energy  of 
si  i  Vdt  u>fi  will  result  in  one  limiting  form  (through  the 
matrix  dement,)  arid  a  vanishing  difference  in  the  free  energy 
leads  tu  the  opposite  limiting  form. 

It  is  also  clear  that  any  r. -1  system,  even  if  modelled 
minimally  with  the  use  of  this  double  harmonic  well  potential, 
must  exhibit  a  complicated  collection  of  changes  in  all  of  the 
quant  it  ie.s  which  we  have  considered.  Thus,  in  fact,  it  snould 
not  generally  be  possible  to  reduce  most  transfers  to  the  limit 
expressed  in  eqn  Nevertheless,  this  limit  is  useful  for 

graphically  showing  how  the  quantity  D  can  assume  values  which 
yield  different  limiting  forms  for  the  incomplete  gamma  functions. 

Ir.  general,  the  q  -  quar.t  i  t  ies  satisfy  the  large-argument 
lm.it  ii,r  the  incomplete  gamma  functions.  We  conclude  this 
section  by  considering  the  two  extreme  limits  D  '  0  and  D  >  3 
tor  the  Lj(.,  matrix  element  as  expressed  by  eqn  (29)-(31). 

First ,  the  large  q-limit  for  the  matrix  element  <bO|V^|bO> 


lim.bOjVhbO.  -  ‘  (ka/kb,<q  +ZbK,]e 


In  the  limit  of  a  symmetric  well  system,  this  simplifies  further, 
k^/k^  »  1.  The  large  D- limit  for  the  two-centre  matrix  element 
is  given  by 


1  tor.  ■  aO  •  Vb  |  b0>  •  -  lb.,  S  ■'  -- i  ♦  k  iV.  1 

D;J  '  I2.r(l+a-/b  * 

(  )\  \ 

+  2  I - i - i  b  j  (k  /k,  )  R.  +  R  ’U'D' 

IHa-Zb-l  '  a  b  ^  “Ij 


+  _ 1 _ n  + 

1 +a - / b • [  l+a‘/b; 


In  contrast,  the  small  D- limit  is  given  by 


-  Z  <k.'kK>'„  +  *h  b 


-  iK-  b-  -  2v0/n. 


For  an  essentially  symmetric  well,  this  reduces  to 


-hwbS00q . 


It  is  worth  noting  in  eqn  t^O)  the  fact  that  although  an 
apparent,  direct  dependence  on  V0  disappears  in  the  limiting  process 
nevertheless,  V0  still  plays  a  role  in  the  matrix  element  through 
the  quantity  D. 

Turning  our  attention  to  the  tunnelling  of  a  small  cation, 
such  as  lithium,  from  one  cage  of  solvent  to  another  in  solution, 
we  find  the  following.  First,  assume  that  the  system  satisfies 
the  limits  for  a  symmetric  well  model.  that  is.  the  initial  and 
final  states  are  physically  the  same  In  this  case,  it  is  easy 
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t  it  find 

LaO.W)  *  -  jr.  »M  '■'q  (*-’) 

If,  for  lithiurn,  we  assume  that  w  *  A30  cm"*  (which  is  a  reasonable 

value  based  on  the  far  infrared  spectrum  of  lithium  in  a  number 

7  °  -1 

of  solvents  )  and  t  -  0.3  A.  then  Lag  »  0.05  kJ  mol  V  As 

pointed  our  in  Part  I*,  this  value  satisfies  the  diabatic  limit 

for  the  transfer.  It  remains  to  be  demonstrated  whether  a 

value  for  the  transfer  matrix  element  Lflg  ^0  which  is  obtained  in 

this  fashion  is  ar.  accurate  estimate.  It  is  clear  that  it  is  a 

reasonable  and  consistent  value  when  it  is  compared  with  the 

values  which  were  generated  in  Part  1  by  means  of  the  use  of  more 

complete  and  accurate  potential  energy  functions. 

Discussion  and  Summary 

The  effects  of  the  environment  on  the  matrix  elements  have 
not  been  specifically  included  here,  as  they  were  in  Part  I.* 

The  reason  is  simple.  As  was  argued  In  Part  I,  the  specific 
consideration  of  the  environmental  contributions  to  the  matrix 
elements  for  the  transfer  could  be  reduced  to  additional  terms 
which  are  of  the  same  general  form  as  which  is  defined  for 
the  reactive  subsystem.  For  the  model  system,  as  we  treat  it  here 
In  terms  of  the  double  harmonic  oscillator  potential,  such 
environmental  contributions  are  implicitly  contained  In  the  values 
of  the  force  constants  which  are  used  for  the  initial  and  final 
states.  Thus,  one  assumes,  the  experimental  values  of  the  force 


.  1C. 

constants,  which  arc  used  in  the  calculation,  contain  the  requisite 
Information  about  the  effect  of  the  environment  n  *  h«  *  .it  rlx  clcrrerir 
for  the  transfer. 

The  experimental  values  of  the  force  constants,  however, 
yield  no  direct  information  which  can  be  used  to  estimate 
the  energy  of  activation  In  order  to  estimate  this  energy,  it  is 
necessary  to  resort  to  the  usual  collection  of  models  of  the 
solvent  for  example,  the  polar  continuum  This  problem  of  the 
calculation  of  the  energy  of  activation  is  important,  it  was 
not  considered  here,  but  shall  be  considered  in  a  separate  paper. 

The  model  which  has  been  developed  in  this  paper  provides 
a  method  which  can  be  used  to  estimate  the  site  of  the  pre- 
exponential  factors  for  rate  constants  for  reactions  which,  one 
presumes,  cake  place  in  the  diabatic  limit  In  order  to  carry 
out  the  calculations,  one  need6  only  the  (observed)  frequencies 
for  the  vibrations  of  the  transferring  species,  polarized  along 
the  transfer-axis,  in  its  initial  and  final  states  and  an  estimate 
of  the  distance  over  which  the  species  must  migrate  It  is  not 
necessary  to  have  complete  knowledge  of  the  specific  solute- 
solvent  interaction  potential  energy  functions. 

As  in  Part  1,  this  work  also  clearly  indicates  that  the 
vibrational  overlap  between  the  states  involved  enters  as  the 
principal  factor  which  determines  whether  a  transfer  is  likely 
or  not.  Indeed,  both  Part  1  and  this  work  suggest  that  a  simple 
estimate  of  the  vibrational  overlap  in  the  ground  state,  eqn  (19). 

/  offers  a  quick  appraisal  of  the  likelihood  that  a  particular 

configuration  (with  its  specific  transfer-distance)  can  contribute 
to  the  rate  in  a  diabatic  limit.  Such  an  estimate  can  be  carried 
,  out  with  the  use  of  the  quantities  a  and  b  which  are  determined 
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with  che  use  of  the  experimental  frequencies. 
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Figure  Caption 

The  model  potential  energy  for  the  asymmetric  double  harmonic 
oscillator. 
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